Spin Dynamics and Resonant Inelastic X-ray Scattering in Chromium with 
Commensurate Spin-Density Wave Order 
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We theoretically investigate spin dynamics and La-edge resonant inelastic X-ray scattering (RIXS) 
of Chromium with commensurate spin-density wave (SDW) order, based on a multi-band Hubbard 
model composed of 3d and 4s orbitals. Obtaining the ground state with the SDW mean-field 
approximation, we calculate the dynamical transverse and longitudinal spin susceptibility by using 
random-phase approximation. We find that a collective spin-wave excitation seen in inelastic neutron 
scattering hardly damps up to ~0.6 eV. Above the energy, the excitation overlaps individual particle- 
hole excitations as expected, leading to broad spectral weight. On the other hand, the collective 
spin-wave excitation in RIXS spectra is almost masked by large spectral weight coming from particle- 
hole excitations with various orbital channels. This is in contrast with inelastic neutron scattering, 
where only selected diagonal orbital channels contribute to the spectral weight. However, it may be 
possible to detect the spin-wave excitation in RIXS experiments in the future if resolution is high 
enough. 

PACS numbers: 75.10.Lp, 75.30.Fv, 75.50.Ee, 78.70.Ck, 78.70.Nx 



I. INTRODUCTION 

After the discovery of iron-pnictide superconductors 
whose parent compounds are antiferromagnetic metal, 
the spin dynamics of itinerant antiferromagnetic systems 
with multiple 3d orbitals has attracted much attention. 
In order to elucidate such spin dynamics, we focus on 
another system, chromium, which is known to show spin- 
density wave (SDW) states and their physical properties 
have been well established theoretically and experimen- 
tally,!-^ 

Chromium and its alloys have a body-centered cubic 
(bcc) structure. There are many first-principles band- 
structure calculations of paramagnetic Cii^^— The Fermi 
surfaces of Cr consist of an octahedral electron pocket 
at the r point (the wave vector ^=(0,0,0)) and an oc- 
tahedral hole pocket at the H point (fe=(27r/a,0,0), be- 
ing a the lattice constant), which arc almost the same 
shape. This leads to a nesting vector that stabilizes a 
SDW order. Exactly speaking, these octahedrons are not 
the same size, and the incommensurate SDW occurs in 
Cr. The magnitude of the nesting vector is slightly less 
than 27r/a (the nesting vector is (27r/a)(0.952,0,0))..^ The 
doping Cr with cither Mn or Fc changes the nesting vec- 
tor into (27r/a,0,0). When the Mn (Fe) concentration x 
in Cri_2;Mn^ {Cvi^^Fe^) exceeds 0.003 (0.02), the com- 
mensurate SDW state is stabilizedi^ 

The spin dynamics of the incommensurate Cr has been 
investigated by inelastic neutron scattering. The in- 
commensurate spin excitations in the low-energy region 
(< 100 meV) have been clarified experimentallj*^ and the- 
oreticallyi^ii^ However, spin dynamics in the high-energy 
region, where the interplay of collective spin-wave exci- 
tation and individual particle-hole excitation is expected 
to occur, has rarely been examined except for the data 



up to 550 mcV with large error bar 

From the theoretical viewpoint, it is necessary to in- 
clude all of 3d orbitals in models for the precise descrip- 
tion of both the collective and individual excitations, as 
was proven in the study of the antiferromagnetic phase of 
iron arsenides^ For the incommensurate SDW state, it 
is not easy to perform band-structure calculations includ- 
ing all of the orbitals, since we need to use a large unit 
cell. In fact, commcnsurability with the vector (27r/a,0,0) 
has been assumed in the previous first-principles calcula- 
tions.— In model calculations for the commensurate 
SDW state, with two-band models used, the detailed 
band structures have been ignoredJ^ Therefore, it is im- 
portant to perform the study of spin dynamics of Cr by 
using a precise band structure, even though the commen- 
surate SDW state is assumed. 

In this study, we investigate Cr with the commensu- 
rate SDW state by using a multi-band Hubbard model 
composed of 3d and 4s orbitals. After a self-consistent 
calculation based on the SDW mean-field approximation, 
we obtain the dynamical spin susceptibility by employing 
random- phase approximation (RPA) . We find a collective 
spin-wave excitation undamped up to ^0.6 eV. Above 
the energy, the excitation overlaps individual particle- 
hole excitations, leading to broad spectral weight. From 
the results of the longitudinal spin susceptibility, we find 
that its spectral weight is mainly distributed above the 
energy of the spin-wave mode. We expect that these 
features may be observed in inelastic neutron scattering 
experiments in the near future. 

Resonant inelastic X-ray scattering (RIXS) tuned for 
the L edge of transition metal has been recognized as a 
powerful tool to investigate not only charge and orbital 
(d-d) excitations but also spin excitation in the energy 
and momentum spacesJ^ By using a fast-collision ap- 
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proximation for the RIXS process, we calculate Cr L-edge 
RIXS intensity and compare the spectrum with inelas- 
tic neutron scattering spectrum. We find large spectral 
weight coming from particle-hole excitations with vari- 
ous orbital channels. This eventually masks the collective 
spin-wave excitation. This is in contrast with cuprates 
where the spin-wave excitation has clearly been observed 
in the Cu i-cdge RlXSJ^ii^ However, since there are 
remnants of the spin-wave excitation, it may be possible 
to observe the excitation in RIXS in the near future. 



II. COMMENSURATE SDW STATE OF 
CHROMIUM 

We consider a multi-band Hubbard Hamiltonian H = 
Ho + Hi^t^ where Hq denotes tight-binding hopping terms 
with five 3d and 4s orbitals and Hint denotes interaction 
terms for the 3d orbitals. Hq is given by 



(1) 



This paper is organized as follows. In Sec. |TI1 the 
multi-band Hubbard model with five Sd orbitals and 4s 
orbital is introduced together with the SDW mean-field 
approximation for the commensurate SDW order param- 
eter. In Sec. mil the dynamical susceptibility is calcu- 
lated within RPA. Predictions to inelastic neutron scat- 
tering are made. In Sec. IIVI RIXS spectra tuned for Cr 
L3 edge are calculated based on a fast-collision approxi- 
mation. The origin of characteristic spectral distribution 
is clarified. Summary is given in Sec. [V] 

I 



where cj ^ ^ creates an electron with orbital ji and spin a 
at site i. The matrix clement t{A.ij] fi, v) represents the 
hopping of an electron between the orbital at site i and 
the V orbital at site j with the distance of j- = r; — rj . 
The values of t{A.ij; fj., v) are obtained by a down-folding 
procedure of the first-principles band structure in the 
paramagnetic phase. We use the first-principles calcula- 
tion package pwSCF (Ref. 20)and WannierQO (Ref. 21)for 
the down-folding procedure. We set the lattice constant 
of the bcc lattice to be a=2.883 Aand take 512 points 
in the first Brillouin zone (BZ) of the bcc lattice. In 
Hq, we take the hopping matrix elements up to the 25-th 
neighbors. 

Hint is given hy21 
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where the indices /i and v run over only 3c? orbitals, and 
U and J are on-site Coulomb and exchange interaction 
for 3d electrons, respectively. 

In order to consider the Hamiltonian in the momentum 
space, we substitute 



c -J-Ve''^''^ 



Ck,, 



(3) 



and its hermitian conjugate into Eqs. ([T]) and ([2|), where 
N is the number of sites. Imposing the commensurate 
SDW order with a ordering vector (5=(27r/a,0,0), we ob- 
tain a mean-field Hamiltonian 
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with the order parameter 



(7) 



The average (■ • • ) is taken at zero temperature in the 
present study. In Eq. constant terms are neglected. 
We can rewrite the momentum fc with feo defined within 
a reduced zone due to the commensurate SDW order: 



k ^ ka + mQ. 
This leads to the follwoing replacement: 

E-E E ■ 

k ko m— 0,1 



(8) 



(9) 



We introduce quasi-particle operators j^{kQ,a) for 
band e, satisfying 

^^i^tJ.,m;e{ko,<j)'^e{ko,(j), (10) 



where ipfj,,m;e{ko,(j) is the eigenvector of i/^F ^^j, 
eigenvalue Eko.a,e- The mean-field Hamiltonian reads 

^""^ = E E ^feo,<..al(fco, Theiko, a). (11) 



We self-consistently solve mean-field equations to de- 
termine the order parameters. From a self-consistent so- 
lution, we obtain Ekg^a-,e and tpp,,m-eikQ,a). In the self- 
consistent calculation, we use 50 x 50 x 50 meshes in the 
first BZ of the bcc lattice. The chemical potential is de- 
termined so that the number of electrons should be 6.0 
per site, i.e., the number of the valence electrons in Cr. 

The values of U and J are unknown. Considering that 
the magnetic moment defined by 



(12) 



■m.fi 



is sensitive to the choice of U and J, we determine their 
value so that the calculated M should be close to the 
observed one (0.6 /iB, being /iB the Bohr magneton) in 
Cr,^ We obtain U=2.5 eV and J=0.1U. Note that there 
are several possible values of U and J, but the conclusions 
presented in this paper are independent of the choice of 
U and J. 

The band structure and the density of states calcu- 
lated by our method are consistent with previous publi- 
cations.—"— We find that the calculated magnetic mo- 
ment, 0.6 fiB^ is composed of 0.48 fiB in orbitals and 
0.12 fiB in eg orbitals. 



III. SPIN AND CHARGE DYNAMICS 
A. Dynamical Susceptibility 

We define the dynamical susceptibility as 

1 X - 



Xl\{q,q',uj) 



N ^ 

k',k" 



dt e 







X ([4',A,a2Wcfc'+q,K,'Ti(i),4"+q',^,^^Cfc",^,^j]), (13) 

where Ck^fi^ait) is the Heisenberg representation of Ck.fi.ai 
and s {s') denotes a spin pair (7iCT2 (c'lCri)- * takes f, i, 
-I- and — corresponding to (t,t), {iA), (lit), and (t,l), 
respectively, and the same for s' . In an abbreviated form, 
q' is omitted for the q = q' case as x('7,w). 

The bare susceptibility is represented with the wave 
function and quasiparticle energies, 

xo' >~^M^q + iQ.^) 

J_ /(-E'po+q,<Ti,e) — /(£'po,o-2,e') 

N ^ 



X ?/'«:,m;e(P0 + O"! )Cm+n;e (Po + 9 + IQ , Cr2) 

X ■4jtrn-e'iP0,<^2)'4't,,7n+n-,e'{PQ,(j[)Sai.a'Ja{.a2, (14) 

where the sum with respect to po runs over the reduced 
zone. We set rj = 0.01 eV. 

We calculate the dynamical susceptibilities within the 
multi-orbital RPA, 
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where the product of susceptibility x and interaction V 
is taken as a matrix product represented in the orbital 
basis such as 

[XoV^x]Af,(9,w) 

E XoxAQ,Q + ^T'Q,^)K[x'Xx'^^.iq + mQ,q,uj). 

(16) 

The nonzero elements of the interaction matrix V are 
shown in Table ID 

The transverse spin susceptibility , the longitudi- 
nal spin susceptibility x^, and the charge susceptibility 
X" are given by 



x+-(g,c^)-)_^x+;(g,c^), (17) 
riq^Lu) = Y.Y. ~ xf.(g,^) j , (18) 
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TABLE I. The non-zero elements of V^f^ , where s = (ai, (72) and s' — (a'l, a'2). Each of the four subscripts takes one of the 3d 



orbitals, 4s not taken into account. 
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and 

x"iq,u^)^Yl E lxlMq,^) + xZ{Q,^)], (19) 
respectively, where s denotes 4- (t) for * =t (i)- 

B. Spin and charge excitations in commensurate 
SDW state 

Figure [TJa) shows the imaginary part of the transverse 
spin susceptibiUty, Imx^ ('Zjw), in the commensurate 
SDW state of Cr from q'=(0,0,0) to (27r,0,0). Here, the 
lattice constant a is taken to be a unit length. We note 
that this q range corresponds to a range from q=27r(- 
1,1,0) to 27r(0,l,0). In Fig. [TJa), there appears disper- 
sion with strong intensity starting from the magnetic 
zone center ^=(271, 0,0). This is a collective spin-wave 
excitation. The collective excitation exhibits a promi- 
nent weight up to around 0.6 eV, but gradually loses 
its sharpness and weight above 0.6 eV. A possible cause 
of this damping is individual particle-hole excitations. 
To investigate this, we plot the particle-hole excitation 
spectrum, i.e., the imaginary part of the bare transverse 
spin susceptibility, Imxo~(QjW), in Fig. [TJc). We find 
that the collective excitation gradually penetrates into 
the particle-hole continuum around 0.6 eV. Therefore, it 
is clear that the damping is caused by the particle-hole 
excitations. A similar origin of damping has been dis- 
cussed for iron arsenides based on the same calculation 
scheme 

The spectral weight distribution and dispersive feature 
in Fig. [He) come from the band structure in the com- 
mensurate SDW phase. Therefore, the threshold energy 
(~0.5 eV) of particle-hole excitation at <7=(0,0,0) may be 
seen through optical process. In fact, the calculated peak 
of the optical conductivity induced by the commensurate 
SDW order— is located at 0.6 eV (not shown here), sim- 
ilar to the threshold energy. This energy is also roughly 
consistent with an observed peak position (^0.45 eV) of 
the optical conductivity for commensurate Cr alloys,— 
although the calculated one shows a slightly larger value. 
The 0.15 eV overestimate of the peak position in our the- 
ory as compared with the experimental value may come 
from neglecting the effect of band renormalization due to 
correlation. This amounts to a factor of 3/4 for the band 
renormalization: the energy scale in our theory should be 



(a) (b) 




FIG. 1. (Color online) Contour plot of the imaginary part of 
spin susceptibility in the commensurate SDW state along the 
(g,0,0) direction, (a) Imx"' (transverse) and (b) Imx^ (longi- 
tudinal), (c) Imx^j (transverse) and (d) Im^o (longitudinal) 
for the bare susceptibility. The magnitude of intensity shown 
by the color bar in (a) is limited up to 20, but the maximum 
intensity is 460 at around q = 2n and w ~ 0.1 eV. We use 
40 X 20 X 20 meshes in the first BZ of the bcc lattice. 



multiplied by 3/4 to make a comparison with experiment. 
The same factor is obtained from the comparison of the 
peak positions of the optical conductivity in the param- 
agnetic phase between theory (1.2 eV)^ and experiment 
(0.9 eV).^ 

The longitudinal spin excitation is shown in Fig. [TJb). 
The strong spectral weight is located above 0.4 eV near 
q=(27r,0,0). We find that the main part of the weight is 
distributed above the collective spin-wave energy shown 
in Fig. [IJa). This means that, if the experimental setup 
of inelastic neutron scattering is properly chosen to de- 
tect both the transverse and longitudinal excitations, a 
broad and less dispersive spectral distribution is expected 
near the magnetic zone center up to 1 eV. By taking 
into account the fact that the magnetic component par- 
allel to the scattering vector does not contribute to the 
scattering, this kind of setup may be achieved for both 
a longitudinal SDW (L-SDW) phase below 122 K and 
a transvers SDW (T-SDW) phase between 122 K and 
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FIG. 2. (Color online) Contour plot of the imaginary part of 
dynamical susceptibility near the BZ center, (a) Transverse 
spin mode, Imx"' , (b) longitudinal spin mode, Imx^, and 
(c) charge mode, Imx" along the (g,0,0) direction. The mag- 
nitude of intensity is shown by the color bar on the right of 
each panel. We use 80 x 20 x 20 meshes in the first BZ of the 
bcc lattice. 



311 K in Crji assuming that the high-energy excitations 
in the incommensurate phases are similar to those in the 
commensurate phase: The momentum transfer q for L- 
SDWshould be set to near (0,27r,0), while q for T-SDW 
should be set to near (27r,0,0) or (0,27r,0). The experi- 
mental confirmation of these features is highly desired. 

The difference of the spectral properties between the 
transverse and longitudinal excitations naturally comes 
from the difference of the bare spin susceptibility as 
shown in Figs.HJc) and (d): The spectral weight of Imxg 
(Fig. [T{d)) near q=(27r,0,0) is distributed higher in en- 
ergy than that of ImxJ~ (Fig. IHc)). 

To see the detailed structures of gapless excitations, 
we investigate the spectra around the T point, where the 
spectral intensity is moderately weak enough to see the 
whole structure clearly. We also note that small q re- 
gion contributes to L-edge RIXS spectra as discussed in 
Sec.[IVl The imaginary part of x"*""; X^j and x" near the 
r point are shown in Fig.[2lja), Fig.[2jb), and Fig.[2ljc), re- 
spectively. The intensity of the transverse mode, Imx"' , 
is larger than that of the longitudinal mode, Imx^, be- 
low 0.5 eV. The velocity of the transverse mode Vs is 
estimated to be 1.8 eVA. This value is consistent with 
a theoretically estimated value based on a band gap (^ 
1.5 eVA)^^ We notice that Vs is smaller than that 
of the charge mode. This is expected from the large 
U limit, where the exchange interaction determining the 
spin- wave velocity, /{U — J), is smaller than the hop- 
ping integral t that determines the energy scale of charge 
motion. 



IV. RESONANT INELASTIC X-RAY 
SCATTERING SPECTRUM 



In this section, we calculate RIXS spectra in Cr tuned 
for the L3 edge. First, we derive the formula of spectral 
intensity within a fast-collision approximation, where the 
dynamical susceptibilities calculated in the former sec- 



tion are included. Next, the calculated RIXS spectra are 
analyzed in terms of spin, charge, and orbital degrees of 
freedom. 



A. Formulation of RIXS spectral intensity 

The L-edge RIXS consists of two processes: the X-ray 
absorption and emission by way of an intermediate state 
with core holes. The X-ray absorption is accompanied 
by the creation of a Cr 2p core hole and an electron in id 
and 4s orbitals. The 2p core-hole state has a hybridized 
spin with the orbital angular momentum due to spin- 
orbit coupling. When the electron relaxed through the 
X-ray emission is not that excited from 2p, the resulting 
orbital occupation of Sd and 4s electrons varies from the 
initial. Note that the spin of the relaxed electron may be 
either up or down, since the core-hole spin is hybridized. 
The final state may thus involve a spin flip. This means 
that we can investigate spin excitations in addition to 
charge and orbital ones through the L-edge RIXS.— 

We introduce a transition operator Dk associated with 
the X-ray absorption and emission processes: 



-\-k.}j..(T 



Pk' 



+ h.c, 



(20) 



where Pk'j.j^ is the annihilation operator of Cr 2p core- 
electrons and cji^~ (e) is the dipole-matrix element, being 
e the unit vector of the polarization of the incoming and 
outgoing X ray. The matrix element is given by 



(21) 



where \2p,j,jz) represents the 2p state with the total 
angular momentum j whose z component is jz , and a) 
represents the 3d and 4s states with orbital /.i and spin a. 
The values of c^'-'^ (e) are estimated by using 2p, 3d and 
4s atomic orbitals. The ratio of the radius part of the 
2p-4s element to that of the 2p-3d element is ^ 0.20, 
implying that the 2p-3d dipole transition is dominating 
in Cr L-edge RIXS. 

By using the second order perturbation involving the 
X-ray absorption and emission processes, the intensity of 
the RIXS spectrum reads 

Imxsiq = fcin - fcout, W = Win " Wout) 



{f\Dl 



1 



^ Win + Ea - Him + iF 



(22) 



where L^o (Ef) is the energy of the initial (final) state 
and 1 /r is the life time of the intermediate state. Him is 
a Hamiltonian that acts on the intermediate state. We 
employ a fast-collision approximation, assuming that 1 /F 
is very short, which corresponds to neglecting the details 
of the energy dispersion in the intermediate state. This 
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FIG. 3. The geometry of scattering process used in the RIXS 
calculation, q — fein — fcout. The scattering plane is a;z-plane, 
and q—{q,0,0). We change q by varying the angle a between 
fcin and /Cout . 



approximation leads to a simple expression of the RIXS 
spectrum given by 

S.S' K^X.ll.U 

X C^fffi (s^out) C^'^J^ (^in)''A'i<?3 (^out)cift7*4 (£in 

(23) 

where £i,i (Sout) denotes the polarization vector of incom- 
ing (outgoing) X-ray. Hereafter, we consider the Cr L3 
edge for resonance, i.e., j = 3/2. 

The geometry for our RIXS calculation is illustrated 
in Fig. [3l The angle between the momenta of incom- 
ing (fcin) and outgoing (fcout) X-ray is a, and the norms 
of these momenta are almost same. We change the mo- 
mentum transfer q — km — fcout by varying the angle 
a. q=(g,0,0) is on the x-axis. The incoming beam is 
either 7r-polarized (parallel to the scattering plane) and 
cr-polarized (perpendicular to the scattering plane). The 
outgoing X-ray is chosen to be a summation of both 
polarizations, taking into account existing experimental 
conditions. The spin is assumed to direct parallel to the 
z-axis, which is perpendicular to the ordering vector Q. 
This geometry corresponds to one of two spin polariza- 
tions in the T-SDW phase between 122 K and 311 K in 
the incommensurate Cr.— We note that, since a charac- 
teristic X-ray energy for Cr L3 absorption is 572 eV, the 
momentum space where Cr La-edge RIXS can access is 
limited to q 0.537r. 



B. Calculated RIXS spectra 

The calculated RIXS spectra for the commensurate 
SDW state are plotted in Fig. |4] for both the tt and a po- 
larizations along the (g,0,0) direction. The high-energy 
part of the spectra above 1 eV is shown in Figs. HJa) 
and (b). The spectra are asymmetric with respect to 




FIG. 4. (Color online) Contour plot of the calculated RIXS 
spectra of Cr for the commensurate SDW state along the 
(g,0,0) direction. High-energy region (1 eV< u <3 eV): (a) 
TT polarization and (b) a polarization. Low-energy region 
(0 eV< u <1 eV): (c) tt polarization and (d) a polarization. 



q=0. This is not caused by the dynamical suscepti- 
bility in Eq. (f^. but is due to the product of dipole 
matrix elements 1^ Decomposing the polarization of the 
outgoing X-ray, we found that the asymmetry is predom- 
inately originated from the 7r-polarizcd outgoing X-ray. 
The change of incident-photon polarization induces the 
shift of the intense spectral cloud at high energy ^ 1 eV 
from (7=— 0.347r for the tt polarization to — O.Stt for the a 
polarization. 

The intensity of the spectra increases with increasing 
energy and shows a broad maximum at a; ~ 2 eV. With 
further increasing energy, the intensity decreases. The 
main contribution to the broad peak comes from inter- 
band excitations between different orbitals as discussed 
below. The broad peak at w ~ 2 eV is expected to appear 
in RIXS experiments for Cr. 

The low-energy region below uj=l eV is shown in 
Figs.m^c) and (d). According to the dynamical spin sus- 
ceptibility shown in Fig. [51 we expect low-energy spin 
and charge excitations staring from uj=0 cV at <j'=(0,0,0). 
However, the present scale of intensity in Figs. |31[c) and 
(d) makes the excitations almost invisible because of the 
presence of high intensity in the high-energy excitations. 
However, the low-energy excitations certainly exist even 
though their intensity is small. 

In order to make clear the presence of the low-energy 
excitations indicated by Fig. [31 we show the line shape 
of the RIXS spectra for both polarizations in Fig. \5\ We 
find dispersive low-energy excitations, starting from near 
w=0 eV and ^=(0,0,0), with the same energy scale as 
the spin and charge modes shown in Fig. ^ although it 
is difficult to identify each mode. 

One may have a question why the spin and charge 
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FIG. 5. (Color online) The line shape of RIXS spectra 
of Cr in the commensurate SDW state for (a) n- and (b) 
CT-polarization. Each line is plotted for momentum transfer 
between (— 2l7i-/35,0,0)-(2l7r/35,0,0) from bottom to top with 
the increments of 7r/35. 



modes are less clear in the RIXS spectra. In order to 
resolve this question, we should notice a difference of the 
suffixes of the dynamical susceptibility between RIXS in 
Eq. (|23p and, for example, the transverse spin mode in 
Eq. p7)) . In Eq. p7)) . only k = \ and ji ^ v, i.e., exci- 
tations within the same orbital, arc taken from the dy- 
namical susceptibility in Eq. while in Eq. k. A, 
/i. and V are independently taken. In order to see the ef- 
fect of the difference, we calculate the RIXS spectra only 
taking k = X and ^ ~ v m. Eq. (|23p and show the results 
in Fig. [6l We find that the spectral weight near 1 eV is 
dramatically suppressed and the low-energy modes be- 
come clear in intensity. This means that the excitations 
between different orbitals dominate the RIXS intensity. 
Therefore, wc may conclude that the collective modes like 
spin-wave excitation cannot be easily observed by L-edgc 
RIXS, if the systems under consideration contain many 
orbital components near the Fermi level. We note that 
this does not happen in the case of single-band Hubbard 
model describing cuprate compounds. 



V. SUMMARY 

We have investigated the spin dynamics and RIXS for 
Cr with the commensurate SDW state by using a multi- 
band Hubbard model composed of Sd and 4s orbitals. 
We have employed a self-consistent calculation of the 
ground state based on the SDW mean-field approxima- 
tion. On top of the ground state, we have calculated 



FIG. 6. (Color online) Contour plot of the RIXS spectra 
obtained by restricting the orbital indices in the dynamical 
susceptibility to only hi = \ and fi = v in Eq. (|23|) . (a) the 
TT polarization and (b) the a polarization along the (q,0,0) 
direction. 



the dynamical spin susceptibility within RPA. By eval- 
uating the hopping integrals from the down-folding pro- 
cedure of first-principles band calculation, the electronic 
states in Cr are properly described by our model. The 
electronic states in the commensurate SDW phase ob- 
tained by the mean-field theory are also consistent with 
those from the previous first-principles band structure 
calculations. Therefore, the calculated dynamical sus- 
ceptibilities reflect realistic particle-hole excitations. It 
is well-known that RPA cannot describe the band renor- 
malization effect due to correlation. The factor of the 
band renormalization is estimated to be roughly 3/4, by 
comparing the peak positions of the optical conductivity 
between theory and experiment as discussed in Sec. lIIIBl 

We have found that a collective spin-wave mode ap- 
pears in the spin-transvers excitation spectrum. The 
collective mode does not damp up to ~0.6 eV. Above 
the energy, the excitation overlaps individual particle- 
hole excitations, leading to broad spectral weight. We 
expect that this feature appears in inelastic neutron scat- 
tering experiments for Cr, even though Cr shows an in- 
commensurate spin excitations below 0.1 eV. When the 
contribution from the longitudinal spin excitation to in- 
elastic neutron scattering increases, the spectral weight 
above the spin-wave mode is expected to be enhanced. 
This may also be detected by future experiments. 

RIXS tuned for the L-edge is a powerful tool to inves- 
tigate not only charge and orbital excitations but also 
spin excitation in the energy and momentum spaces. By 
using a fast-collision approximation for the RIXS pro- 
cess, we have calculated Cr La-edge RIXS intensity. We 
have found an asymmetric spectral distribution along the 
(g,0,0) and (— g,0,0) directions. This asymmetry comes 
from the dipole matrix element between 2p-core orbitals 
and 3d orbitals. We have also found large spectral weight 
coming from intcrband excitations between different or- 
bitals. This eventually masks collective spin-wave exci- 
tations within the same orbital. Even though the exci- 
tations are weak in intensity, it may be possible to de- 
tect them if the experimental resolution of RIXS is high 
enough. 
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